Origin of the roughness exponent in elastic strings at the depinning threshold 
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Within a recently developed framework of dynamical Monte Carlo algorithms, we compute the 
roughness exponent £ of driven elastic strings at the depinning threshold in 1 + 1 dimensions for 
different functional forms of the (short-range) elastic energy. A purely harmonic elastic energy leads 
to an unphysical value for £. We include supplementary terms in the elastic energy of at least quartic 
order in the local extension. We then find a roughness exponent of £ ~ 0.63, which coincides with the 
one obtained for different cellular automaton models of directed percolation depinning. The quartic 
term translates into a nonlinear piece which changes the roughness exponent in the corresponding 
continuum equation of motion. We discuss the implications of our analysis for higher-dimensional 
elastic manifolds in disordered media. 



The competition between elasticity and disorder is a 
central theme of current research in statistical physics. 
The pinning of flux lines in a type-II superconductor |Q , 
the motion of a charge density wave 0, an interface in 
magnets and several realizations of surface growth || are 
all governed by these two antagonistic mechanisms, one 
trying to smooth the surface, the other striving to distort 
it. These two mechanisms are already at work in the sim- 
plest such system, a one-dimensional elastic string in a 
two-dimensional medium, which today is far from being 
solved. In this paper, we report progress in our under- 
standing of the zero-temperature motion of this system 
at the depinning threshold, defined by the critical driv- 
ing force f c . Above f c , the elastic string flows with finite 
velocity, while it is pinned for forces / < f c . We consider 
the problem on the lattice, but keep full contact with the 
continuum description. 

Previously [0, we showed that the dynamical Monte 
Carlo method can be reconciled with the continuum 
equation of motion approach if extended, non-local, 
moves are allowed for. In fact, earlier Monte Carlo work 
|^],[| had been hindered by pathologies of the commonly 
used local move set which leads to an infinite critical force 
for unbounded disorder. Many workers in the field 
had circumvented these problems by rather considering 
cellular automaton models. These models are very useful. 
However, it is normally impossible to identify the differ- 
ential equation which results in the continuum limit. In 
contrast, the Monte Carlo dynamics is derived from an 
energy function, and it satisfies detailed balance. We de- 
veloped a method which finds with great ease the critical 
string, i. e. the string which is pinned at the critical force 
f c . In Q), we showed that the critical force and the crit- 
ical string are completely independent of all the details 
of the (non-local) Monte Carlo algorithm. Both from a 
conceptual and a practical point of view, the situation is 
thus much better controlled. In this paper, we are con- 



cerned with the statistical properties of critical strings 
for different functional forms of the (short-range) elastic 
energy. 

Specifically, we consider an elastic string h — 
{/i'}i=o,... ,l moving on a finite lattice of size L x M with 
periodic boundary conditions in both directions. For con- 
creteness, we take our random potential V(i,j) at each 
site to be made up of uncorrelated Gaussian variables 
with unit variance. The energy of a string h l in presence 
of an external driving force / is given by 

L 

E (tf) = £ {V(i, hi) - fh\ + E d (Al)} . (I) 

i=i 

Here, the (short-range) elastic energy i? e j is a function of 
the local extension A' = h\ +1 — h\. Periodic continuation 
is implied, such that, for large driving forces /, the string 
h* keeps winding around the finite lattice. The Variant 
Monte Carlo method presented in |i) allows to compute 
the critical string h c for an arbitrary local convex elas- 
ticity. 

The string's roughness exponent £, in the thermody- 
namic limit, is defined by 

(h$-hff~c<mBk\i-j\*. (2) 

Following pc[ , we obtain £ by computing as a function 
of system size L the mean square elongation 

W 2 (L) := < {h c - < h c >) 2 >. (3) 

In eq.(||), < h c >= j- J2i hi: while the overbar stands for 
an average over the disorder. 

As a first step, we show in figure ^ the mean square 
elongation W 2 {L) as a function of L for a harmonic elas- 
tic energy i? e j(A) = A 2 for system sizes (L x M) up to 
(1000 2 ). The data are very well fitted by a straight line of 
slope 2.34, which indicates that the roughness exponent 
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is £ ~ 1.17. It is now well understood that a line with 
C > 1 cannot represent a physical string |ll). In fact, 
the thermodynamic limit breaks down in this case, and 
the structure of the elastic string is described by a size- 
dependent constant in eq.(||) \h\ — h^\ 2 ~ A 2 (L)\i — j\ 2 , 

with a diverging mean square extension A 2 (L) in the 
limit L^oo ||. 
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FIG. 1. Mean square elongation W 2 (L) as a function of 
system size L for a quadratic elastic energy I? e j(A) = A 2 (up- 
per curve) and for the model with metric constraint |A| < 1 
and a = (c/eq.@). The interpolating lines correspond to 
roughness exponents of £ = 1.17 and 0.63, respectively. 

The same behavior has been observed in the contin- 
uum limit (i — > a;) of the lattice model eq. ([!]), where a 
harmonic elastic energy yields a term ~ V 2 in the cor- 
responding time evolution equation (quenched Edwards- 
Wilkinson equation): 

= f + r 1 {x,h) + 2a 1 V 2 h. (4) 

at 

Here, the random force r](x, h) is the negative derivative 
of the random potential V in eq.([D). The e q.(^) has been 
much studied by analytical methods 13|-15| as well as 
by direct numerical simulation |[l~0)J]~T| , |l6| , |l 7| and the ex- 
istence of a roughness exponent in excess of one is now 
well accepted. From the numerical work, we expect a 
value of £ ~ 1.15 while Chauve et al Jl5| obtained 
£ ~ 1.2 ± 0.2 from a two-loop functional renormaliza- 
tion group calculation. Our numerical data in figure ^ 
thus establish perfect agreement of the lattice with the 
continuum framework. 

The unphysical roughness of the elastic string de- 
scribed by the quenched Edwards- Wilkinson equation 
has led many authors |l(],[ll|] to conclude that a 
1— dimensional string necessarily develops overhangs and 



islands, which go beyond a description by a single- valued 
function h(x,t). Others PJl8| have attempted to intro- 
duce and to justify nonlinear terms in the time evolution 
equation (||) (cf below). 

In this paper, we study elastic strings which can be 
described by a function h(x,t). We pursue a approach 
based on the analysis of the energy eq.(|l|) : We argue that 
a diverging extension A is contradictory with the replace- 
ment of £Li(A) by the A 2 term, dominant only at small 
A|. Correspondingly, the quenched Edwards- Wilkinson 
equation stems from a small gradient expansion |3| , and is 
incompatible with this situation. We find that £ changes 
and becomes physical (and universal) if terms beyond 
lowest order in A are kept. A single higher-order con- 
tribution to the elastic energy is found to be important. 
It generates a nonlinear piece in the corresponding con- 
tinuum equation of motion, whose origin has not been 
identified before. 

We first consider an energy function, where the har- 
monic potential is cut off by a metric constraint: 



^el( A ) 




(5) 



By construction, the constraint in eq.(^|) forces the rough- 
ness exponent to be £ < 1. Our data in the lower part of 
figure |l] evidence however a much more interesting fact, 
namely an exponent £ ~ 0.63. 

This exponent coincides with the one obtained in many 
numerical simulations of directed percolation depinning, 
as obtained with a variety of cellular automata models. 
Very interestingly, both the cellular automaton (model 
B) of Sneppen [[7| and the rule proposed by Tang and 
Leschhorn ||] explicitly provide for a local constraint 
which limits the value of |A|, as we can do in eq.(||). 
Notice that a whole family of Monte Carlo rules corre- 
spond to each energy function and that the critical 
strings coincide for all its members. 

As mentioned in the beginning, we strive to keep full 
contact between the Monte Carlo dynamics on the lat- 
tice and the continuum description. This can be achieved 
more easily if we consider smooth (but non-harmonic) 
functions E^A) of which eq.(||) should come out as a 
limiting case. To do so, we note that the elastic energy 
must be a symmetric function of the local extension A. 
We make the Ansatz to write it as a power series in A 2 



£ el (A) = aiA 2 + a 2 A 4 



(6) 



where we assume all coefficients a, to be positive. In 
this way we assure convexity of the function £J„i(A), 
which is a necessary prerequisite for our numerical algo- 
rithm. By direct computation of the roughness exponent 
in a variety of cases, we have found the same exponent 
£ ~ 0.63 whenever a term of at least quartic order in A 
was present. Consider, in figure 0, our data for different 
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choices of the coefficient a\ and a 2 in eq.(||). As can be 
seen in the figure, the roughness exponent is in all cases 
Q ~ 0.63, a value which is not changed by powers A 6 
and higher. This value of £ is physically acceptable, as it 
yields a finite elastic energy per link of the critical string. 
We conclude that the elastic energy E e \(A) contains only 
relevant terms of order A 2 and A 4 . This insight allows 
us to directly consider the continuum limit. 
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FIG. 2. Mean square elongation W 2 (L) as a function of 
system size L for elastic energies incorporating different quar- 
tic terms. The interpolating lines correspond to roughness 
exponents C, = 0.63, as for the metrically constrained curve in 
figure jjj. 

In this limit, the quenched Edwards- Wilkinson equa- 
tion is modified by the negative gradient of the quartic 
term 



d_ 
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h(x, t)=f + rj(x, h) + 2oi V 2 /i + 12a 2 V 2 /i(— ) 2 . 
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Our two observations: 



• including a term A 4 yields a physically satisfactory 
roughness exponent £ < 1; 

• powers A 6 and higher do not influence the value of 

C; 

can now be interpreted in the light of the standard hydro- 
dynamic scaling argument (cf ||). This argument sup- 
poses the string to be self-affine, so that a scale transfor- 
mation x — > bx corresponds to a transformation h — * b^h. 
In the hydrodynamic limit (6 — > 00), higher derivatives 
(corresponding to higher-order terms in the elastic en- 
ergy eq.^)) are negligible compared to lower ones if the 
roughness exponent satisfies £ < 1. Notice that the 



hydrodynamic scaling argument is self-consistent: The 
roughness exponent has to satisfy £ < 1 even after drop- 
ping higher order terms. The naive application of the 
argument would lead us to neglect also the last term 
in eq.(^|). This would bring us back to the quenched 
Edwards- Wilkinson equation, which describes a string 
which is not self-affine. 

The discrepancy between the values of the roughness 
exponents C — 1-17 and £ ~ 0.63 has usually been inter- 
preted as due to the presence of two distinct universality 
classes for interface growth in disordered environments 
P|, one belonging to the quenched Edwards- Wilkinson 
equation (eq.([|)) Qjl^-KL the other represented by the 
cellular automata models (7|-|9| as well as by experiments 
on directed percolation depinning ||. We would argue 
that £ = 1.17 is a specific result, valid for the harmonic 
potential only, and strongly dependent on the behavior 
E el (A) for |A| -► 00. 

Numerical work by Amaral et al fisj had detected the 
presence of a nonlinear relevant term in cellular automata 
that was initially believed to be caused by the A(^) 2 
term in the quenched KPZ equation p0| 

8 8h 

-h(x, t) = f + V (x, h) + v^h + A/2( — ) 2 . (8) 

However, in the KPZ equation this term is of kinematic 
origin and vanishes at the depinning threshold |l^,^l|,^2[ , 
where the critical string is a purely static object. It was 
later believed [|l8|,^l) that anisotropies in the disorder dis- 
tribution might be responsible for such a nonlinear term. 

In our system, an anisotropy is clearly absent. Rather, 
we have shown here that a roughness exponent £ ~ 0.63 
is naturally generated either by a higher order term in 
the elastic energy (c/cq.(^J)) or, equivalently, by a metric 
constraint as in eq.(|g). Such a constraint is also present 
in cellular automata models. Without it, we generate an 
unphysical elastic string with £ ~ 1.17. The nonlinear 
term we introduced into eq.(Q) is the most relevant one 
which can be generated at the microscopic level within 
the static description of an energy function eq. (|]) . Note 
that such a description is sufficient below the depinning 
threshold, where viscous or inertial effects are naturally 
absent. We suspect that the nonlinear term may gener- 
ate other, more relevant, terms upon coarse graining (for 
related work, cf ^|). This point will certainly have to 
be studied analytically, and goes beyond the scope of this 
paper. 

The quenched Edwards- Wilkinson equation has also 
been studied for the driven motion of d-dimensional in- 
terfaces in a d + 1 dimensional target space. The d = 4 
dimensional exact, stable, solution has been used as a 
starting point for perturbative renormalization group cal- 
culations in dimension d = 4 — e p^-|l5||. The fact that, 
in d = 1, the roughness exponent £ > 1 seems to indicate 
that, by continuity, our nonlinear term will be present 
even in higher dimensions. There are indications of at 
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least logarithmic divergencies even in d = 2 (fij. Clearly, 
much effort will still be needed in order to completely 
understand this system. 

In conclusion, we have studied in this paper the ef- 
fect of the elastic energy on the statistical properties of a 
string at the depinning threshold. For a harmonic energy 
E c y ~ A 2 , we find a roughness exponent £ familiar from 
numerical and analytical work on the quenched Edwards- 
Wilkinson equation. This exponent indicates that the 
string is unphysical, as the mean square extension of the 
string A 2 diverges in the thermodynamic limit. We then 
studied different elastic energies, either by cutting off the 
harmonic function with a metric constraint on |A| or by 
using a more general (differentiable) Ansatz for £^,i(A). 
In both cases we regularized the string and obtained a 
physically sound roughness exponent C — 0.63. The low- 
est order addition to the elastic energy, of quartic order 
in A, generates a nonlinear piece in the corresponding 
continuum time evolution equation. 

Although we do not know at present how this new non- 
linear term (which preserves te symmetry h — > — h) will 
behave under coarse graining, we can be quite sure that 
it is correct on a microscopic level and gives a consis- 
tent thermodynamic description of a driven elastic string 
at the depinning threshold. It will generate a roughness 
exponent of £ ~ 0.63 in the continuum. 

Finally, we note that the Monte Carlo methods intro- 
duced in Q] have allowed to obtain the critical string even 
for large systems, and for general convex energy func- 
tions. This yields the completely transparent transition 
between the two situations with £ = 1.17 and £ ~ 0.63, 
of which only the latter is universal and has physical rel- 
evance 
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